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Abstract. We present a comprehensive study of the symmetries of the generating 
functionals of generic Langevin processes with multiplicative colored noise. We treat 
both Martin-Siggia-Rose-Janssen-deDominicis and supersymmetric formalisms. We 
summarize the relations between observables that they imply including fluctuation re- 
lations, fluctuation-dissipation theorems, and Schwinger-Dyson equations. Newtonian 
dynamics and their invariances follow in the vanishing friction limit. 
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1. Introduction 

The stochastic evolution of a classical system coupled to a quite generic environment 
can be described with the Langevin formalism [TJ [2] and its generating functional, the 
Martin-Siggia-Rose-Janssen-deDominicis (MSRJD) path-integral [21 IH IS]- In many 
cases of practical interest the effect of the environment is captured by an additive 
white noise and its memory-less friction, Brownian motion being the paradigmatic 
example [TJ. Nevertheless, there are many other interesting instances in which the 
noise is multiplicative and colored, and the friction effect is consistently described 
by a memory kernel coupled to a non-linear function of the state variable. Such 
Langevin equations appear in many different branches of physics (as well as chemistry 
and other sciences). In magnetism, the motion of the classical magnetic moments of 
small particles is phenomenologically described by the Landau-Lifshitz-Gilbert equation 
in which the fluctuations of the magnetic field are coupled multiplicatively to the 
magnetic moment [7j. Many other examples pertain to soft condensed matter; two 
of these are confined diffusion, in which the diffusion coefficient of the particle depends 
on the position via hydrodynamic interactions [8], and the stochastic partial differential 
equation that rules the time-evolution of the density of an ensemble of N Brownian 
particles [Sj. In a cosmological framework, these equations arise as effective equations 
of motion for the out of (although close to) equilibrium evolution of self-interacting 
quantum fields in which the short-wave length modes serve as thermal baths for longer 
wave-length modes with slower dynamics pH]. Such type of fluctuations may yield a 
priori unexpected results such as noise induced phase transitions in systems in which 
the associated deterministic potential does not exhibit any symmetry breaking |llj . 

In order to better understand these processes it is useful to distinguish cases in 
which sources of fluctuations and dissipation can be different. On the one hand, the 
noise and friction term can have an 'internal' origin, such as in diffusion problems. On 
the other hand, the stochastic fluctuations can be due to an 'external' source [12J. In 
the former cases one usually assumes that the variables generating the noise and friction 
are in equilibrium and the terms in the Langevin equation associated to them are linked 
by a fluctuation-dissipation theorem. In the absence of non-conservative external forces 
the Boltzmann measure of the system of interest is a steady state of its dynamics. In 
the latter cases noise and dissipation are not forced to satisfy any equilibrium condition 
and this translates into the possibility of having any kind of noise and friction terms. 
For concreteness we shall focus on the first type of problems and only mention a few 
results concerning the latter. 

In treatments of the examples mentioned in the first paragraph, the delicate double 
limit of vanishing fast variables relaxation time and noise correlation time is often taken. 



These lead to a first order stochastic differential equation with multiplicative white noise. 
Its interpretation in the Ito, Stratonovich or other sense requires a very careful analysis 
of the order of limits, see e.g. [13] and references therein. In the main part of this paper 
we shall keep both time-scales finite and thus avoid the subtleties encountered in the 
double vanishing limit. 

In this manuscript, we identify a number of symmetries of the MSRJD generating 
functional of inertial Langevin processes with multiplicative colored noise. One 
symmetry is only valid in equilibrium. The corresponding Ward-Takahashi identities 
between the correlation functions of the field theory lead to various equilibrium relations 
such as stationarity, fluctuation-dissipation theorems [IJJ[T5] or Onsager relations. Away 
from equilibrium, the symmetry is broken giving rise to various out of equilibrium 
fluctuation relations [16]- [28]. Another symmetry holds for generic out of equilibrium 
set-ups and implies dynamic equations coupling correlations and linear responses. It 
allows to express, in particular, the linear response in terms of correlations without 
applying a perturbing field [30]- [JT]. 

We are aware of the fact that some of the results in this manuscript - especially, 
in the limit of additive noise - were already known and we do our best to attribute 
them to the authors of the original papers or review articles. Still, the presentation that 
we gradually develop in this article allows one to go beyond the simple cases and treat 
the multiplicative non-Markov processes with the same level of difficulty Moreover, we 
discuss in greater detail than previously done the transformation of the measure and 
several Jacobians, and the domain of integration of the fields in the path-integral. The 
importance of dealing with colored noise, and to treat the transformation of the fields 
in the complex plane, is enhanced by our purpose to extend this analysis to quantum 
dissipative problems. These results will be presented in a separate publication [42J. 

The organization of the paper is the following. In Sect. [2] we give a short review 
of Langevin equations with additive and multiplicative noise. Section [3] presents the 
MSRJD functional representation of Langevin equations. In Sect. H]we deal with the 
equilibrium symmetries while in Sect. |5]we treat the out of equilibrium ones. In both 
Sections we discuss supersymmetric formulations. We conclude in Sect. [6j 

2. The Langevin equation 

We consider a 0-dimensional field ip ( e -9- a particle at position ip) with mass m driven 
by a force F and in contact with a thermal bath in equilibrium at inverse temperature 
p. The initial time, to, is the instant at which the particle is set in contact with the 
bath and the stochastic dynamics 'starts'. We call it to = ~T and without loss of 
generality we work within a symmetric time- interval t G [— T, T] . The extension to 
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higher dimensional cases is straightforward. Our conventions are given in |Appendix A 



2.1. Additive noise 

The Langevin equation with additive noise is given by 

EQ{[iP],t)=mi>{t)-F{[4>],t)+ I duj{t,u)iP{u) = £{t), (1) 

J-T 

with ip(t) = dip{t)/dt and ip{t) = d 2 t(j(t)/dt 2 . The force can be decomposed into 
conservative and non-conservative parts: F([^],t) = — V'(ip(t), \(t)) +f([ip],t). V is a 
local potential the time-dependence of which is controlled externally through a protocol 
A(t). V denotes the partial derivative of V with respect to ip. f([ip],t) collects all the 
non-conservative forces that are externally applied, f ([?/>], t) is assumed to be causal in 
the sense that it does not depend on the future states of the system, tp(t') with t' > t. 
Furthermore, we suppose that f ([?/>],£) does not involve second - nor higher - order 
time-derivatives of the field ip{t). The last term in the left-hand-side (lhs) and the 
right-hand-side (rhs) of the equation model the interaction with the bath. These two 
heuristic terms can be derived using a model [I3j 01] in which the bath consists in a set 
of non-interacting harmonic oscillators of coordinates qi that are bilinearly coupled to 
the state variable of the system of interest ip. The function 7 is the retarded friction 
[7(i, £') = for t' > t] and the noise £ is a random force taken to be a Gaussian process. 
This assumption is quite reasonable, for instance, for a Brownian particle much bigger 
than the surrounding particles of the bath, its motion being the result of a large number 
of successive collisions, which is a condition for the central limit theorem to apply. Since 
we assume the environment to be in equilibrium, 7(£,£') is a function of t — t' and the 
bath obeys the fluctuation-dissipation theorem of the 'second kind' [45J: 

<£(*)>* = 0, <£(t)«0>€ = /r 1 r(* - , (2) 

where (...)§ denotes the average over the noise history. We introduced the symmetric 
kernel T(t — t') = j(t — t') + 7(t' — t) = T(t' — t). The white noise limit, in which the 
bath has no memory, is achieved by setting / ~f(t — t') = 7o<5(£ — £') with 70 > the friction 
coefficient. The Langevin equation then takes the more familiar form 

eq(m,*) = mm - F{M,t) +70 m = m ■ (3) 

Newtonian dynamics, for which the system is not in contact with a thermal bath, are 
recovered by simply taking j(t) = V(t) = at all t. Out of equilibrium environments 
can be taken into account by relaxing the condition between the noise statistics and the 
friction kernel T(t - t') = j(t - if) + j(t' - t). 



2.2. Multiplicative noise 

We generalize our discussion to the multiplicative noise case in which the Gaussian noise 
£ is coupled to a state-dependent function M'(ip). The Langevin equation reads 



EQ(ty>],£) = mi/>{t) - F{[ip],t) + M'(xp(t)) / du-f(t-u)M'(ij(u))ij(u) 

=M'm))m- (4) 

This equation can also be shown by using the oscillator model for the bath and a 
nonlinear coupling of the form M(ijj) ^V c^ where q are coefficients that depend on 
the details of the coupling and M(ip) is a smooth function of the state variable with 
M (0) = 0. By a suitable renormalization of 7, one can always achieve M' (0) = 1. For 
reasons that will soon become clear, we need to assume that M'{ip) ^0 V^. These 
assumptions can be realized with functions of the type M(ip) = ip + L(ip) where L is a 
smooth and increasing function satisfying L(0) = 1/(0) = 0. The complicated structure 
of the friction term takes its rationale from the fluctuation-dissipation theorem of the 
second kind that expresses the equilibrium condition of the bath. This equation models 
situations in which the friction between the system and its bath is state-dependent. £ 
has the same statistics as in the additive case, see eq. (J2J). The Langevin equation for 
the additive noise problem is recovered by taking M(ip) = if). 

2.3. Initial conditions 

The Langevin equation is a second order differential equation that needs two initial 
values, say the field and its derivative at time — T. We shall use initial conditions drawn 
from an initial probability distribution P; (ip(—T),ip(—T) ) and average over them. The 
initial conditions are not correlated with the thermal noise £. In the particular case in 
which the system is prepared in an equilibrium state, p is given by the Boltzmann 
measure. 

2.4- Markov limit 

Langevin equations are often given in the Markov limit in which they appear to be first 
order stochastic differential equations. Second and higher order time-derivatives as well 
as non-local terms such as memory kernels are not allowed. In other words, the effect 
of inertia is neglected (Smoluchowski limit) and the bath is taken to be white. This 
is justified in situations in which the two associated time scales are sufficiently small 
compared to all other time scales involved. Concretely, the resulting equation is derived 
by using an adiabatic elimination procedure that consists in integrating over the fast 
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variables of the system (the velocities) and of the bath. However, this double limiting 
procedure requires a careful analysis. 

The physics of the resulting equation may depend on how the relaxation time 
associated to inertia compares with the correlation time of the noise before sending 
the two of them to zero. In cases in which the latter is much larger than the former, 
the limiting stochastic equation should be interpreted in the sense of Stratonovich |47j . 
The RHS of eq. 03]) is given a meaning by stating that ip in M'(if)(t)) is evaluated 
at half the sum of its values before and after the kick. Conversely, when the inertia 
relaxation time is much larger than the noise correlation time, the limiting equation 
should be interpreted in the Ito sense [35]. In this scenario, the rule is that M'(ip(t)) 
is evaluated just before the kick £(£). When the noise is additive the two conventions 



are equivalent (see |Appendix B.2[ ) for all practical purposes. However, they are not 
for processes with multiplicative noise [12]. For these it is possible to rewrite the Ito 
stochastic equation in terms of a Stratonovich stochastic equation by adding an adequate 
drift term to the deterministic force - and be allowed to use the rules of conventional 
calculus. The Fokker-Planck equation associated to the Markov process does not depend 
on the scenario and the Boltzmann distribution is a steady state independently of the 
convention used. However, the action of the generating functional acquires extra terms 
depending on the discretization prescription [8] 132] • 

In this article, we decide not to cope with the Markov limit and, unless otherwise 
stated, we keep the inertia of the system in our equations (m ^ 0) and we use a colored 
noise with a finite relaxation time. 

3. The generating functional 

The generating functionals associated to the equations of motion (JTJ) and ((3]) are given 
by the Martin-Siggia-Rose- Janssen-deDominicis (MSRJD) path-integral. In this Section 
we recall its construction for additive noise [3] and we extend it to multiplicative noise 



by using a continuous time formalism. In |Appendix B| and |Appendix C| we develop a 



careful construction in the discretized formulation. 



3.1. Action in the additive noise case 



The Langevin equation (JT]) is a second order differential equation with source £. The 
knowledge of the history of the field £ and the initial conditions ip(—T) and ip(—T) is 
sufficient to construct ip(t). Therefore, the probability P[tp] of a given ip history between 
— T and T is linked to the probability of the noise history -P n [£] through 

PWpty] = P n [£] V[Z] Pi (V(-T), ^(-T)) dV(-T) (ty(-r) (5) 



implying 

P[^] = P n [Ec#]] \Jty\\ Pj (>(-T),^(-T)^ 
where J"^] is the Jacobian which reads, up to some constant factor, 



J[ip] = det, 






det, 



£Eq(M,m) 
fo/>(v) 



Jo 



(6) 



(7) 



det [...] stands for the functional determinant. We introduced the notation Jb|/0] f° r 
future convenience and we shall discuss it in Sect. 13.31 After a Hubbard- Stratonovich 
transformation that introduces the auxiliary real field i[>, the Gaussian probability for a 
given noise history to occur reads 



P^]=AT- 1 V$]e 



- fdu ii/>(u)£(u)+| ffdudv i-tp(u)f3~ 1 r(u-v)iip(v) 



(8) 



with the boundary conditions ip(—T) = ip(T) = and where all the integrals over time 
run from —T to T. In the following, unless otherwise stated, we shall simply denote 
them J . M is a infinite constant prefactor that we absorb in a re-definition of the 
measure T>[ip]. Back in eq. (E]) one has 



with the MSRJD action functional 

S[iP,$] = lnpf^(-T),^(-T)) - fduij>(u)EQ([iP],u) 

dudv iij){u) (3~ l T(u - v) vij>(v) + In \J [ip]\ 
The latter is the sum of a deterministic, a dissipative and a Jacobian term, 



(9) 



1 

+ 2 



(10) 



with 



S det [ijj^] = InPi U(_r),V(-T) 



du i^(u) mip(u) — F([i})],u_ 
S diss [ip, ip] = / d« ty(u) /dv t(w - v) [/3 _1 i^(w) - ^(v) 



(11) 

(12) 

S" 1 * 5 * takes into account inertia and the forces exerted on the field, as well as the measure 
of the initial condition. S dlss has its origin in the coupling to the dissipative bath. In 
the white noise limit, j(t — t') = 7o<5(£ — if), the dissipative action naively simplifies to 
S dmB [ip, ip] — 70 / du iip{u) /3 _1 r0(u) — ift(u) (see Sect. 12 .41 for additional details on this 
limit). 
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Integrating away the auxiliary field i\) yields the Onsager-Machlup action 
functional [50]. However, we prefer to work with the action written in terms of ip 
and k/> as this is the form that arises as the classical limit of the Schwinger-Keldysh 
action used to treat interacting out of equilibrium quantum systems [HJ |5T], that we 
shall analyze along the same lines in |42] . 



3.2. Action in the multiplicative noise case 

To shorten expressions, we adopt a notation in which the arguments of the fields and 
functions appear as subindices, ip u = ip(u), 7 U _„ = 7 (it — v), and so on and so forth, 
and the integrals over time expressed as J" = j_ T du . 

In the case of the Langevin equation (J3J) with multiplicative noise, the relation (JBJ) 
is modified and reads 

'EqUb 



P[m = Pn 

with the Jacobian 

J[ip] = det 



_M'(ip) 

8Eq u [^}/M'(iIj u 
8ip v 



\J[i;}\ Pi^-T^-r) 



det, 



S u - 



and the generalization of the definition of j7o in eq. ((7j) to the multiplicative case: 

~<JEq u [V>] M"(ip u ) 



Jo 



det, 



to th 

Eq u hp}5 u - V 



Jo 



(13) 



(14) 



(15) 



8ip v M'(ip u ) 

The construction of the MSRJD action follows the same steps as in the additive noise 
case, complemented by a further transformation of the field vip (->■ iijj M'(ip), the Jacobian 
of which cancels the first determinant factor in the RHS of eq. (TT4|) . Therefore, the 
MSRJD action reads 



S[^ } iP] = lnPi(^_ T ,^_ T ) - / i^E Qu [^] 

J u 

+ \ f JifaM'fa) P- X T U „ V M'(4> v )ii> v + In IJoMI , (16) 

with Jq defined in eq. (ITBI) . The deterministic part of the action is unchanged compared 
to the additive noise case and the dissipative part is now 



S dlss [V,V] = tyu M'^ u ) lu _ v M'^ v ) 



13 hj) v - ip v 



(17) 



3.3. The Jacobian 



In |Appendix C| we prove that the Jacobian j7o is a field-independent positive constant 
for Langevin equations with inertia and multiplicative colored noise. One can therefore 
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safely include the Jacobian contribution in the normalization. However, we decide to 
keep track of this term and represent it as a Gaussian integral over Grassmann conjugate 
fields c and c*. 



J [iP] = Jv[c,c* 



e S J {c,c*,i,] 



with 



and the boundary conditions: c(—T) = c(—T) = c*(T) = c*(T) = 0. Plugging in the 
Langevin equation (jlj), we arrive at 



^'^//<^-/<^ E ^< < 19 ' 



S J [c,c*,if>] = I I c;. 

' U J V 




SF lib] 
md 2 J u -v ~ -^f 1 + M'^ u )d ulu . v M'{^ v 

Oip v 

"1 c: ^w Wm " FuM1cu ' (20) 

The Grassmann fields c and c* that enter the integral representation of the determinant 
are known as Faddeev-Popov ghosts and can be interpreted as spinless fermions. The 
two-time fermionic Green function defined as 

(c$cv)*7= [vic,c*}c* t c t ,e sJ ^^ , (21) 



is related, by use of Wick's theorem, to the inverse operator of — fe' — M , 5tV Eq [ip t ] 8 t _ t i . 
{clcti)$j inherits the causality structure of the latter and it vanishes at equal times as 
long as the Markov limit is not taken {i.e. all fermionic tadpole contributions can- 
cel): (c^Ct')sJ = for t > t'. The last statement can be easily verified by considering 
the discretized version of S^ (see Appendix B.3.3 and Appendix C) and by checking 



that the diagonal terms of the inverse operator vanish in the continuous limit. S^ 
only involves combinations of the form c*c, i.e. it conserves the fermionic charge and 
( c t)s J = ( c t)s J = 0- This implies that S J [c,c*,ip] and more generally the MSRJD 
generating functional (at zero sources) are invariant under the following field transfor- 
mation 

Tj(a) = \ °l * aC _\\ \/aeC*. (22) 

| c* h-> a c* , 

The Jacobian of the transformation is trivially equal to one and the measure T>[c, c*} is 
left unchanged. One has 7j((x)Tj(j3) = Tj(aj3). 

The total MSRJD action given in eq. f fl6|) can be written equivalently as a functional 
of ip, ip, c and c* provided that the path-integral measure is extended to the newly 
introduced fermionic fields: 

Sty, 4>, c, c*\ = S det ty, $} + S diss ty, $\ + S J [c, c\ V>] . (23) 
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S.Jf.. Observables 

3.4-1- Measure. We denote (...) the average over the thermal noise and the initial 
conditions. Within the MSRJD formalism, the average is evaluated with respect to the 
action functional S[ip, ip] or S[ip, if), c, c*] and we use the notation ( ... ) 5 : 

(...) s = [v[if),i>] ...e 5 ^ = [v[i>J,c,c*] ...e s ^ c ^ . (24) 



3-4-2. Local observable. The value of a generic local observable A at time t is a function 
of the field and its time-derivatives evaluated at time t, i.e. a functional of the field if) 
around t, A([ip],t). Unless otherwise specified we assume it does not depend explicitly 
on time and denote it A[^(£)]. Its mean value is 

(Am)]) = (mt)})s . (25) 

3.4-3. Time-reversal. Since it will be used in the rest of this work, we introduce the 
time-reversed field ip by ip(t) = if)(—t) for all t. The time-reversed observable is defined 

as 

A r (M,t) = A([$,-t). (26) 

It has the effect of changing the sign of all odd time-derivatives in the expression of 
local observables, e.g. if A[^(i)] = d t ip(t) then A r [-i/>(£)] = —d t ip(t). As an example for 
non-local observables, the time-reversed Langevin equation ([T|) reads 

EQ r ([^],t) = m^(t)-F r ([^],t)- / dwy(u-t)i>(u). (27) 

Notice the change of sign in front of the friction term that is no longer dissipative in 
this new equation. 

3-4-4- Two-time correlation. We define the two-time self correlation function as 

c{t,t') = m)^')) = m)^'))s. (28) 

Given two local observables A and B, we similarly introduce the two-time generic 
correlation as 

C {AB} {t,t') = (Am)\Bm')])s, (29) 

The curly brackets are here to stress the symmetry that underlies this definition: 

C{AB}(t,t') = C{BA}(t',t). 
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3.4-5. Linear response. If we slightly modify the potential according to V(ifj) (->■ 
V(ip) — fif,ip, the self linear response at time t to an infinitesimal perturbation linearly 
coupled to the field at a previous time t' is defined as 



R{t,t') = 



S{m) 



W) 



Hm)s lu] 



(30) 



U=o 



It is clear from causality that if t' is later than t, (ip{t))s[u] cannot depend on the 
perturbation U(t') so R(t,t') = for t' > t. At equal times, the linear response R(t,t) 
also vanishes as long as inertia is not neglected (m ^ 0]j. More generally, the linear 
response of A at time t to an infinitesimal perturbation linearly applied to B at time 
if < tis 



RAB{t,t ) = 

with V(V>) >->■ V(V>) - / B 5 



(5<A[^(t)]> 



*/e(tO 



*(^^W]>5[/ £ 



/s=0 



*/fl(*') 



(31) 



/b=0 



where B is local. 



5. 5. Classical Kubo formula 

By computing explicitly the functional derivative 5/8 f^ in the path integral generating 
functional, we deduce 



H--)s[u] 



su(t) 



(•• 



5 S[ip,ip, 0,0*;^] 



U=o 



5 Uit) 






(...ms + (...^M^ mt ) )s 



(32) 



M'm)) 

The first term in the RHS comes from the functional derivative of S det . The second term 
comes from the Jacobian term expressed with the fermionic ghosts, S^ , and vanishes 
identically (see the discussion on the equal-time fermionic Green function in Sect. I3.3J) . 
One has 

6 ( 1 )s[U] 



mt))± 



suit) 







U=o 



l${t)ty{i!)) s 



S 2 (l)s iU 







(33) 



(34) 



U=o 



SUit) 5U(t>) 

From the definition of the linear response, eq. (130~I) . we deduce the 'classical Kubo 
formula' |15] 

R(t,t') = (mW)) s . (35) 

| In the double limit of a white noise and m — > 0, the equal-time response can slightly violate the 
causality principle depending on the order in which the limits are taken. In the Ito scenario it vanishes 
whereas in the Stratonovich one it has a finite value. 



14 

The linear response is here written within the MSRJD formalism as a correlation 
computed with an unperturbed action. The causality of the response is not explicit, 
nevertheless following the lines in Ref. [8] one can check it is built-ir|§|. Because of this 
expression, the auxiliary field ip is often called the response field. Observe that we have 
not specified the nature of the initial probability distribution Pi nor the driving forces; 
eq. ( )35l) holds even out of equilibrium. 

Similarly, by plugging eq. (125)) into eq. ( 13"TJ) . we obtain the classical Kubo formula 
for generic local observables: 



R AB (t,t') = (A[tfj(t)] 



SS[ip,ip,c,c*;fB] 



Sf B (t>) 



)s 

/s=0 



wwl /d»i*»)*3ga> s 



^«,;.,^ dB[Uf)\. 



n=0 z t \ i 

This formula is valid in and out of equilibrium and allows us to write the response 
functions associated to generic observables (e.g. functions of the position, velocity, 
acceleration, kinetic energy, etc.) as correlators of ip, ip and their time derivatives. For 
example if B is just a function of the field (and not of its time-derivatives), only the 
n = 0-term subsists in the above sum, yielding 

R AB (t,f) = ( A imm')^^)s . (s?) 

As another example, if one is interested in the response of the acceleration A[-0(t)] = 
dfip(t) to a perturbation of the kinetic energy B[ip(t)] = ^m(d t ip(t)) 2 one should 
compute 

R AB (t,t') = m(dfy(t)d t ,ii>(t')d t ^(t')) s ■ (38) 

Furthermore, it is straightforward to see that within the MSRJD formalism we can 
extend all the previous definitions and formulas to A being a local functional of the 
auxiliary field: A[^(t)]. For example, if A[^(t)] = iip(t) and B[ip(t)] = ip(t), we obtain 
the mixed response 

^(t,t') = (if(tMt')) s = 0, (39) 

where we used eq. ([3T 



In general, a multi-time correlator involving iip(ti) vanishes if ii is the largest time involved. 
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4. Equilibrium 

In this Section we focus on situations in which the system is in equilibrium. We 
identify a field transformation that leaves the MSRJD generating functional (evaluated 
at zero sources) invariant. The corresponding Ward-Takahashi identities between 
the expectation values of different observables imply a number of model independent 
equilibrium properties including stationarity, Onsager relations and the fluctuation- 
dissipation theorem (FDT). These proofs are straightforward in the generating 
functional formalism, demonstrating its advantage with respect to the Fokker-Planck 
or master equation ones, when the environment acts multiplicatively and has a non- 
vanishing correlation time. We shall report soon [42J on the extension to the quantum 
case where the Keldysh action also exhibits a non-trivial symmetry for equilibrium 
dynamics. Similarly to the classical case, this symmetry leads to the quantum FDT. 

4-1. The action 

Equilibrium dynamics are guaranteed provided that, apart from its interactions with 
the bath, the system is prepared and subjected to the same time-independent and 
conservative forces (F = —V). In such situations, the initial state is taken from the 
Boltzmann probability distribution 

InPi^-r^-r) = -0WfoM - ^Z , (40) 

where 'H[ipt\ = ^ m W + ^C0t) is the internal energy of the system, and Z is the partition 
function. The Langevin evolution of the system in contact with the bath can be put in 
the form 

l)C[ ' '" J + M'(il> t ) I lt - u M'{ip u )i> u = M'{ip t )& , (41) 



with C[if) u ] = |w0y — V(ip u ) being the Lagrangian of the system. In this equilibrium 
set-up, the deterministic part of the MSRJD action functional reads 



:»U*r-/. ./.I _ oojU. 1 u. cr / / :„?. 6C i^ 



S aet [ijj, </>] = -PHM-t] -hiZ+ i^ u 




tyu 



= -13 {hmlt T + V{1>-. T )\ - hxZ - J'4 U [mi> u + V'tyj] . (42) 

The dissipative part of the MSRJD action functional remains the same, see eq. ( TT71) . 
As discussed in Sect. 13. 3[ the Jacobian Jo enters the action through the constant term 
In j7o or it can be expressed in terms of a Gaussian integral over the ghosts fields c and 
c*. In this case, its contribution to the action reads 



3 J [c, c*, V] = / / c* [md 2 u 5 u . v + M'(^ u )9 u7lt _M'(^)] c v 






M'(V„ 



M'(^) 
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(43) 



^.S. Symmetry of the MSRJD generating functional 



We shall prove that JP^, ■0, 



c, c 



*1 pSI^^v'] 



is invariant under the equilibrium field 



transformation: 



T 

'cq 



^ u H> V-« , C u !->• C*_ u 

iip u H- i^_ u + (3d u tp- u , c* (-)■ -c- 



(44) 



This transformation is involutary, % q T cq = 1, when applied to the fields ip or hp 
and the composite field c*c. It does not involve the kernel 7 and includes a time- 
reversal. It is interesting to reckon that the invariance is achieved independently by the 
deterministic (S det ), the dissipative (S dlss ) and the Jacobian (S J ) contributions to the 
action. This means that it is still valid in the Newtonian limit (7 = 0). The detailed 
proof that we develop here consists of two parts: we first show that the Jacobian of 
the transformation is unity, then that the integration domain of the transformed fields 
is unchanged. Afterwards we show that the action functional S[ip,i/),c,c*] is invariant 
under T eq . 



4-2.1. Invariance of the measure. The equilibrium transformation T eq acts separately 
on the fields ip and iip on the one hand, and the fields c and c* on the other. The 
Jacobian J~ eq thus factorizes into a bosonic part and a fermionic part. The bosonic part 
is the determinant of a triangular matrix: 



J ( 



cq 



det 



KM) 



5(T eq ip, T eq ip) 



det 



1 

uv 



5V- 



Sip— 1 





84>- 



(det J- [Su+v])' 



and it is thus identical to one [53]. It is easy to verify that the fermionic part J7" c 
as well. 



f 

cq 



4-2.2. Invariance of the integration domain. Before and after the transformation, the 
functional integration on the field if) is performed for values of ip t on the real axis. 
However, the new domain of integration for the field ip is complex. For a given time t, 
ipt is now integrated over the complex line with a constant imaginary part —i(3d t ipf One 
can return to an integration over the real axis by closing the contour at both infinities. 
Indeed, the integrand, e 5 , goes to zero sufficiently fast at ip t — > ±00 for neglecting the 
vertical ends of the contour thanks to the term /3 _1 7o(i?/>f) 2 in the action. Furthermore, 
the new field is also integrated with the boundary conditions ip(—T) = ip(T) = 0. 

The equilibrium transformation leaves the measure V[c,c*] unchanged together 
with the set of boundary conditions c(—T) = c(—T) = c*(T) = c*(T) = 0. 
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4-2.3. Invariance of the action functional. The MSRJD action functional 

S[ip,ip,c,c*} = S [i/),i/j] + S dlss [ip,ip} + S^ {c, c*,ip) is invariant term by term. The 
deterministic contribution given in eq. ( 142]) satisfies 

S^T^Xj] = IuP^tJt) - j[4_ u + f3d u ^ u }[md 2 u ^ u + V'(iP- u )] 

J u 

= lnPi^T, fa) ~ / i^K + V'Oj u )} +p ^p u [m^ u + V'(^ u )] 



= lnPi^T, ifr) - / iiK + H^)] + /3 / 9„7#„] 

J u J u 

= S dct {ipJ], (45) 

where we used the initial equilibrium measure hip (?/>,■?/>) = —/37i[il>] —\nZ. In the first 
line we just applied the transformation, in the second line we made the substitution 
u I—)- —u, in the third line we wrote the last integrand as a total derivative the integral 
of which cancels the first term and creates a new initial measure. 

Secondly, we show that the dissipative contribution 5' dlss [^,^], defined in eq. (JT2J1 . 
is also invariant under the equilibrium transformation. We have 

J U J V 

= j[4 u -^ u \ [M'(^ u ) lv - u M'(ij v )f3-Hip v 

Ju J V 

= S dta [V,$. (46) 

In the first line we just applied the transformation, in the second line we made the 
substitution u \-¥ —u and in the last step we exchanged u and v. 

Finally, we show that the Jacobian term in the action is invariant once it is expressed 
in terms of a Gaussian integral over conjugate Grassmann fields (c and c*). We start 
from eq. (143]) 

*J [ /eqC, /eqC , 'eq^J 

c_ M [md 2 J u . v + M'{il>- u )d u j v - v M'(ilf- v )] c% 




+ / C- u 



-V"(^ u ) + 



M "^-M^ + 




M'(^_ u ) ur - M'(^_ u ] 
c* v [md 2 J v . u - M'{ip u )d ulv - u M\ip v )] c u 

M"(ip, 



-V'^-u) 



-v"M + 



M'^ u ) ^ + 



M'(^y 



-vm 



S J [c,c*,^} 



(47) 



In the first line we just applied the transformation, in the second line we exchanged the 
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anti-commuting Grassmann variables and made the substitutions u i— > —u and v i— > —v, 
in the last step we used d v j v - u = —d v ^ u - v and exchanged u and v. 

4-3. Ward-Takahashi identities 

We just proved that equilibrium dynamics manifest themselves as a symmetry of the 
MSRJD action and more generally at the level of the generating functional. This 
symmetry has direct consequences at the level of correlation functions. If A is a generic 
functional of ip and ip, it implies the following Ward-Takahashi identity 

(A[M] ... ) s = (A[%^,%^} ... ) s . (48) 

This identity leads to all the possible equilibrium relations between observables as we 
shall now describe in the following. These relations can be proven without using the 
MSRJD path integral formalism, however our point is to show that the symmetry is 
able to generate all the equilibrium relations without using any other ingredient. 

4-4- Stationarity 

In equilibrium, one expects noise-averaged observables to be independent of the time 
to at which the system was prepared (in our case to = ~~ T). One-time dependent 
noise- averaged observables are expected to be constant, (A^]) = ct, and two-time 
correlations to be time-translational invariant: (A^^-B^/]) = ft-t 1 - Similarly, one 
argues that multi-time correlations can only depend upon all possible independent 
time-differences between the times involved. These statements have been proven for 
additive white noise processes using the Fokker-Planck [54] formalism. The use of the 
transformation % q allows one to show these properties very easily for generic Langevin 
processes. 

One-time observables. Taking A = 1 and letting B be a generic local observable, 
the equal-time linear response vanishes, i?As(t,t) = 0. Using the classical Kubo 
formula ( 136% 

RAB(t,t) = (f^d^ t ^-} s = 0. (49) 

?1=0 4 ^* 

Applying the transformation % q , we find 

RabM ^ <£«-. ||^) S + ?<£*»*-. f^)s • (50) 

n=0 * r l n=0 * r l 

The LHS and the first term in the RHS vanish identically at all times. One is left with 
the second term in the RHS that simply reads (d t B r [ip_ t ]) = d t (B T [ip_ t ]) = 0, proving 
that all one-time local observables are constant in time. 
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Two-time observables. Because we just showed that (A[?lj(t)]) is constant in 
equilibrium, the response Rab^, t'), see its formal definition in eq. ( 13TJ) . can only be a 
function of the time- difference between the observation time and the time at which 
the perturbation is applied. Therefore it can be written in the form RAB{t,t') = 
f(t—t')9(t—t f ). We shall see in Sect. I4.7l that the fluctuation-dissipation theorem relates, 
in equilibrium, the linear response RAB(t,t') to the two-time correlation C{AB}(t,t') 
implying that this last quantity is also time-translational invariant. 

Similarly, (n + l)-time correlators can be proven to be functions of n independent 
time-differences because they are related, in equilibrium, to responses of n-time 
correlators that are time-translational invariant. 

4-5. Equipartition theorem 

Let us consider the local observables A[-0(t)] = d t ip(t) and B[ip(t)] = ip(t). The linear 
response is RAB(t,t') = (d t ip t iif;t>) s — d t {4>ti4>t') s an d we recognize d t R(t,t'). Using the 
field transformation 7~ eq , we find 

d t R(t,t') = dt<jl)-ti1>-t>)s + f3{d t ip-tdt4-t')s 

= dt(il>-tv4>-f)s + (3{d t *Ptdt'*Pt'}s ■ (51) 

If t > t', the first term in the RHS vanishes by causality. Considering moreover the 
limit t! — y t~ the lhs is 1/m as we shall show in Sect. 15.2.31 Finally, we recover the 
equipartition theorem for the kinetic energy 

(3m((d t ^ t f) = 1 . (52) 

4-6. Reciprocity relations 

If we use % q in the expression (1291) of generic two-time correlation functions, we have 

(A[^ t ]B[^])s = (A[^ t ]B r [^ t/ ]) s , (53) 

reading 

C {A B}(t,t / ) = C {Al B r }(^- t ')- ( 54 ) 

In cases in which A and B have a definite parity under time-reversal: 

C{ab}(t) = C{ab}(|t|) if A and B have the same parity, 
C{ab}{t) = -C {A b}{-t) otherwise. 
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4-7. Fluctuation- dissipation theorem (FDT) 

4-7.1. Self FDT. Applying the transformation to the expression ( 135)) of the self 
response R(t, t') we find 

(4>d4>ti) s = (% q tpt%q^t') s = (4>- t i4>-t>)s + P(ip-td t >^-t>)s , (55) 

and we read 

R(t, t') = R(-t, -t') + f3d t >C(-t, -t') (56) 

that, using the equilibrium time-translational invariance, becomes 

R(t) - R(-t) = - /3d T C(-r) , (57) 

where we set r = t — t'. Since C{r) is symmetric in r by definition, this expression can 
be rewritten, once multiplied by 0(r), as 

R{t) = - e{T){3d T C(r) . (58) 

Equation (|58|) is the well-known fluctuation-dissipation theorem. It allows one to predict 
the slightly out of equilibrium behavior of a system - such as the irreversible dissipation 
of energy into heat - from its reversible fluctuations in equilibrium. 

4-7.2. Generic two-time FDTs. We generalize the previous FDT relation to the case 
of generic local observables A and B. Applying the equilibrium transformation 7^ q to 
expression ( 1361) of the linear response RAB(t,t') 



n=0 4 r n=0 ' rt 



+^wE <rv-, fi^ 



9 3pVtf 

n=0 l 



5 



n=0 d W 

Applying once again the transformation to the last term in the RHS yields 

(^]f>» f^)s = (Mi>-t\Y. d ^-t' d -f^)s + Pd t ,(A[i, t ]B[i, t ,})sXm 

which reads 

Rab(t) - R A M-r) = ~ Pd T C {AB} (r) . (60) 

By multiplying both sides by 0(t) we obtain the FDT for any local A and B 

Rab(t) = - &(r)/3d T C {A B}(r) - (61) 
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4-8. Higher- order FDTs: e.g. 3-time observables 

We give a derivation, via the symmetry of the MSRJD formalism, of relations shown 
and discussed in, e.g. j5l], within the Fokker-Planck formalism for stochastic processes 
with white noise. 

4-8.1. Response of a two-time correlation. We first look at the response of a two-time 
correlator to a linear perturbation applied at time t\ 



R(t3,h',ti) = 



(62) 



In the MSRJD formalism, it can be expressed as the 3-time correlator 

R{hM\h) = tytAtjiti) s ■ (63) 

Causality ensures that the response vanishes if the perturbation is posterior to the 
observation times: R(t 3 ,t 2 ;t 1 ) = if ti > max(t 2 ,t 3 ). We assume without loss of 
generality that t 2 < t 3 . Under equilibrium conditions, the response transforms under 

/eq as 

R(ts,t 2 ;ti) = (ip-t^^ip^s + f3d tl (^-t 3 i/j- t2 ip- tl }s ■ (64) 

Multiplying both sides by Q(t 3 — ii) and transforming once again the last term in the 
RHS, the last equation can be written in the form 

{Pd tl (ifj t3 tfj t2 tfj tl )s if h < t 2 < t 3 , 
R{-t 3 , -t 2 ; -h) + (3d tl (ij t3 ij t2 ij tl ) s if t 2 <h<t 3 , (65) 
if t 2 < t 3 < h . 

4-8.2. Second order response. Let us now look at the response to a perturbation at 
time ti of the linear response R(t 3 ,t 2 ): 

In the MSRJD formalism, it can be expressed as the 3-time correlator 

R{h\h,tx) = {ttatytaty^S ■ (67) 



R(t3'i t 2 ,ti) = 



(66) 



It is clear from causality that the response vanishes if the observation time is before the 
two perturbations: R(t 3 ;t 2 ,ti) — if £3 < min(ti, t 2 ). The response transforms under 



/eq &S 



R(t 3 ; t 2 , h) = R(-t 3 ; -t 2 , -h) + f3d tl R(-t 3 , -t i; -t 2 ) 

+ 0dt 2 R(-t 3 , -t 2 ; -h) + pdtM^-t^-t^-t^s ■ (68) 
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Let us assume without loss of generality that t\ < t%. Using causality arguments and 
applying once more the equilibrium transformation to the remaining terms, 

if t 3 < h < t 2 , 

/?(;/,: /,./,.)= { +(3d tl R{t 3 ,t 1 ;t 2 ) iitx<t 3 <t 2 , (69) 

(3d tl R(t 3 ,t i; t 2 ) ift 1 <t 2 <t 3 . 

4-9. Onsager reciprocal relations 

Rewriting twice eq. ( 160]) as 

Rab(t) - Raa(-t) = ~ Pd T C {AB} (r) , (70) 

R BA (-r) - R BtAt (t) = f3d T C {BA} (-r) = f3d T C {AB} {r) , (71) 
and summing up these two equations with r > 

R AB (r) = R BiAi (t) . (72) 

These equilibrium relations, known as the Onsager reciprocal relations, express the 
fact that the linear response of an observable A to a perturbation coupled to another 
observable B can be deduced by the response of B r to a perturbation coupled to A r . 

4-10. Supersymmetric formalism 

4-10.1. Generating functional. The generating functional of stochastic equations with 
conservative forces admits a supersymmetric formulation. This has been derived and 
discussed for additive noise in a number of publications [SSJ EU \57\ . We extend it here 
to multiplicative non-Markov Langevin processes (see [58] for a study of the massless 
and white noise limits). To this end, let us introduce 9 and 8*, two anticommuting 
Grassmann coordinates, and the superfield 

*(*, e, n = m + c *(t) e + e* C (t) + e*e (${ t ) + c *(t) C (t) m;(^)) ) • 

The MSRJD action S [see eq. ( 123|) ] has a compact representation in terms of this 
superfield: 



with 



" '-'susy ' ^susy 1 \' "-*) 

SSy[*] = -I 3 [ d9 d9 * 0*6V,[^{-T, 6, 6*)] - In Z + fdt dd d9* £[*] , (74) 
S^[^} = \ ffdT'dT M(*(T')) D (2) (T', T) M(*(T)) , (75) 
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%\ty\ = ^m^ 2 + V(^) and C[^f] = \mi> 2 — V{^>). In the second equation above we used 
the notation T = (t, 8, 9*) and dT = dt d6 dd*. The 'dissipative' differential operator is 
defined as 

D< 2 >(T', T) = j(t' - t)5(6*' - 9*W - 6) (2/r 1 ^: + sl|, |) , (76) 

where sig^ is a short notation for 26-^ — 1. It is equal to 1 if there is a 6 factor in the 
right and to -1 otherwise. D^ 2 ^ can be written as 

D ( 2 ) ( T , ; T ) = 7 ^/ _ t ^ e *i _ 0*^0/ _ ^ (p B _ B ^ ^ ( 77 j 

with the (covarianttjjl) derivatives acting on the superspace: 



D = l- D = ^-»i- (78) 

that obe}0 {D, D} = -^ and {D, D} = {D,D} = 0. In the white noise limit the 
dissipative part of the action simplifies to 

O*] = \ / dT M (*( T )) D(2) (T) M(*(T)) , (79) 

with the 'dissipative' differential operator 

D (2) (T) = 70 fa- 1 -^ + slg e |) = 7o (DD - DD) . (80) 

This formulation is only suitable in situations in which the applied forces are 
conservative. The Jacobian term S J contributes to both the deterministic (S^ y ) and 
the dissipative part (S^y) of the action. 

4-10.2. Symmetries. In terms of the superfield, the transformation Tj(at) defined in 
eq. f )22|) acts as 

Tj(a) = V(t, 6, 6*) h> #(*, a- l 9, a6*) Va G C* , (81) 

and leaves the action S[^f], see eq. (175)) . invariant. The transformation % q given in 
eq. ( )44l) acts as 

7; q = *(*, 0, 0*) m- *(-i - £0*0, -r , 0) , (82) 

and leaves the action S[^f], see eq. (173)) . invariant. 

The action S[&] given in (173)) has an additional supersymmetry generated by 

* s w< ^s^i- (83 » 

|| Covariant in the sense that the derivative of a supersymmetric expression is still supersymmetric. 
% Therefore the "J 2 term in £[*] can be written in terms of covariant derivatives as ({D, D}^) . 
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that obey {Q, Q} = f and {Q,Q} = {Q, Q} = {D, Q} = {D, Q} = {D, Q} = 

{D, Q} = 0. Both operators Q and Q are thus nilpotent and {Q, Q} is the generator 
of the Lie sub-group. They act on the superfield as 

e e * Q ^ = * + e*Q^ , e eQ ^ = * + eQ^ , (84) 

where e and e* are two extra independents Grassmann constants and 

Q* = c + ^ + c*cf^), (85) 



QH> = - /TV - 9* \B-Hil> - d t ip + p-^cj^j - 9*9 d t c* . (86) 

Expressed in terms of superfield transformations, S[^f] is invariant under both 

V(t,9,9*)^V(t,9,9* + e*) (87) 

and 

V(t, 6, 9*) i->. *(t + e9\ 9 + p~ x e, 9*) . (88) 

Here again, the invariance of the action is achieved independently by the 
deterministic (S det ) and the dissipative (S dlss ) contributions. We would like to stress 
the fact that the presence of the boundary term accounting for the initial equilibrium 
measure of the field ip as well as the boundary conditions for the fields hp, c and c* are 
necessary to obtain a full invariance of the action. 

4-10.3. BRS symmetry. The symmetry generated by Q is the BRS symmetry that 
generically arises when a system has dynamical constraints (here we impose the system 
to obey the Langevin equation of motion). Applying the corresponding superfield 
transformation in (^f(t,9,9*))s gives 

(*(*, 9, 9*)) s = (*(£, 9, 9*) + e*Q*(t, 9, 9*))s , (89) 

and therefore (Q9(t,6,6*)) 8 = 0. This leads to 

(*) s =o, ^ +c ^^ui )s=0 - (90) 

Applying the transformation inside the two-point correlator ( 1 J/ (t , 9 , 9*)^ (t' , 9' , 9*')) $ , 
we find (QV(t,0,0*)V(i;,0',0*'))s + (t,G,0*) <B- (t',9',9*') = 0. This leads to identify 
the two-time fermionic correlator as being the (bosonic) linear response: 



R(t,t') = ^ t 



hp t > +c*,Ct'- 



)s = (4c t ) s . (91) 



M"(W) 

Corroborating the discussion in Sect. I3.3[ this tells us that (c*c t >)s (and more generally 
the fermionic Green function (c*c t ')s^) vanishes for t > t' and also for t = t' provided 

+ e and e* are independent of the coordinates 8 and 0* . 
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that the Markov limit is not taken. Using this result, the second relation in f[9"Uj) now 
yields (i*0 t )s = 0. 

4-10.4- FDT. The use of the symmetry generated by Q on (\l/(£, #, #*))s gives, 

(c* t ) s = , (ty t - /3«9 4 ^>s = . (92) 

By use of (iipt)s — (which was a consequence of the BRS symmetry), the 
second relation becomes d t (ipt)s — 0. This expresses the stationarity and can be 
easily generalized to more complicated one-time observables, A(if)), by use of the 
supersymmetry in {A(^)) S - 

The use of the symmetry generated by Q on a two-point correlator of the superfield 
reads 

(*(t, 9, 6*)V{t', 9', 6*'))s = (* (t + e9*,9 + f3e, 6*)V{t' + e9*\ 9' + 0e, 6*')) s , 

implying, amongst other relations, 



(A 



„ , , M"(if> t ) 

iibt' — pot'ibt 1 + c t ct———— 



<c r )5 = . (93) 



As discussed in Sect. I4.10.3[ (c$cti)gj- vanishes for t > t' . Therefore, the term in c\c t 
disappears from eq. ( 1931) and the FDT is obtained by multiplying both sides of the 
equation by ©(£ — t') 

R(t, t') = pdt?C(t, t')Q{t - t') . (94) 

4-11. Link between T cq and the supersymmetries 

It is interesting to remark that both supersymmetries (the one generated by Q and the 
one generated by Q) are needed to derive equilibrium relations such as stationarity 
or the FDT. All the Ward-Takahashi identities generated by the combined use of 
these supersymmetries can be generated by % q but the inverse is not true. The 
supersymmetries do not yield relations in which a time-reversal appears explicitly such 
as the Onsager reciprocal relations. 

It is clear from its expression in terms of the superfield, eq. (!82j) . that the 
equilibrium transformation % q cannot be written using the generator of a continuous 
supersymmetry. However, the transformation % q can be formally written in terms of 
the supersymmetry generators as 

7; q = ^^nHe^^, (95) 

where II is the time- reversal operator (t i— )■ —t), S exchanges the extra Grassmann 
coordinates {6 \- > —9* and 9* \-¥ 9) and the generator Q is defined in terms of Q and Q 
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as 

Q=-fi9*9{Q,Q} = -fi9*e^. (96) 

5. Out of equilibrium 

We now turn to more generic situations in which the system does no longer evolve in 
equilibrium. This means that it can now be prepared with an arbitrary distribution and 
it can evolve with time-dependent and non-conservative forces f . 

We first show that the way in which the symmetry T eq is broken gives a number of 
so-called transient^ fluctuations relations [H)]-[2HJ ES]- Although fluctuation theorems 
in cases with additive colored noise were studied in several publications [23j-[26j, we are 
not aware of similar studies in cases with multiplicative noise. 

We then exhibit another symmetry of the MSRJD generating functional, valid 
in and out of equilibrium. This new symmetry implies out of equilibrium relations 
between correlations and responses and generalizes the formula? in [30]- [37] obtained for 
additive white noise. Finally, we come back to the equilibrium case to combine the two 
symmetries and deduce other equilibrium relations. 

5.1. N on- equilibrium fluctuation relations 

5.1.1. Work fluctuation theorems. Let us assume that the system is initially prepared 
in thermal equilibrium with respect to the potential V(ip, \-tW The expression for the 
deterministic part of the MSRJD action functional [see eq. (ITT]) ] is 

S det [^, fa A, f] = - m{[4>-T], A_r) - In Z(A_ T ) 



iVv 



m^ u + V'^ w \ u )-{ u [fa , (97) 



where T-t([ip t ],K) = \ m '*P'i + V(ipt,^t)- The external work done on the system along 
a given trajectory between times — T and T is the sum of the work induced by the 
non-conservative forces and the one performed through the external protocol A: 

W[V;A,f] = U u f u [V] + [d u \ u d x V(i> u ,\ u ) . (98) 

J u J u 

The equilibrium transformation % q does not leave S det invariant but yields 

S det [</>, fa A, f] h* S det [V, fa A, f r ] + /3 A T T - f3W[ij- A, f r ] , (99) 

* As opposed to steady-state fluctuation relations the validity of which is only asymptotic, in the limit 
of long averaging times. 

t( This is in fact a restriction on the initial velocities, ip—T, that are to be taken from the Boltzmann 
distribution with temperature (3 , independently of the positions ip-T- The distribution of the latter 
can be tailored at will through the A dependence of V. 
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or equivalently 

S det [</>, ^; A,f] + /3AJF - /W[</>; A,f] k> S dct [</>, ^; A,f r ] . (100) 

S det [ijj^; A, f r ] corresponds to the MSRJD action of the system that is prepared (in 
equilibrium) and evolves under the time-reversed protocol X(u) = A (—it) and external 
forces f r ([?/>],ii) = f([V>], —u). AF T is the change in free energy associated to this time- 
reversed protocol: (3AF T = — \nZ(X(T)) + In Z(X(—T)) = — (3AJ 7 between the initial 
and the final 'virtual' equilibrium states. The dissipative part of the action, S dlss , 
is still invariant under 7~ eq . This means that, contrary to the external forces F, the 
interaction with the bath is time-reversal invariant: the friction is still dissipative after 
the transformation. This immediately yields 

e^(A[^,^]e-^[^ A - f ]) 5[A , f] = (A[T cq ^r eq ^) s[ - XM (101) 

for any functional A of ip and if). In particular, for a local functional of the field, -A[^(t)], 
it leads to the relation ED 



e /3A.F 



(^(t)]e-^ A > f ]) 5[A , f] = (AM-t)])s { xM , (102) 



or also 



r~ T (A[il>(t)]B[il>(l!)]e-Wto*A) 



= (AM-t)]Brm-t')]) s[m . (103) 

Setting A[ip, ijj\ = 1, we obtain the Jarzynski equality [18] 

^(e-^sM = 1 • (104) 

Setting A[%l),ijj\ — S(W — W[ip; X, f]) we deduce the Crooks fluctuation theorem [20l[T9] 

P(W) = P r (-W) e^ w - Ar) , (105) 

where P(W) is the probability for the external work done between — T and T to be 
W given the protocol X(t) and the non-conservative force f ([?/>],£). P T (W) is the same 
probability, given the time- reversed protocol A and time- reversed force f r . The previous 
Jarzynski equality is the integral version of this theorem. 

5.1.2. Fluctuation theorem. Let us now relax the condition that the system is 
prepared in thermal equilibrium and allow for any initial distribution P x . We recall 
the corresponding deterministic part of the MSRJD action functional given in Sect. [31 
eq. (HU) 

S dct [^J]= lnPi(>(-T),^(-T)) 

- fdu v4){u) \mtj){u) + V'{tp{u), X(u)) - i([if>],u)] . (106) 
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The transformation % q does not leave S dot invariant but one has 




S d >,^;A,f]-S^S dct [^;A,f r ], 


(107) 


with the stochastic entropy S = — 


hiP^iT),-^)) -InP^i-TliPi-T)) 


-PQ. 



The first term is the Shannon entropy whereas the second term is the exchange entropy 
defined through the heat transfer Q = AH - W[i/j;\,i]. AH = H([ip{T)}, A(T)) - 
T-L([ip(— T)], A(— T)) is the change of internal energy. The dissipative part of the action, 
S dlss , is still invariant under T cq . This immediately yields 

(A[^]e- 5 > s[A , f] = (A[%^,%J\) sl U) ( 108 ) 

for any functional A of ip and ip. Setting A[i/), ip] = 1, we obtain the integral fluctuation 
theorem (sometimes referred as the Kawasaki identity) 

1 = ^W^f] > (109) 

which using the Jensen inequality gives (S)s\ib$-\t] — 0> expressing the second law of 
thermodynamics. Setting A[if), ip\ = 5(( — S) we derive the fluctuation theorem [171 [19] 

P(C) = P r (-C)e c , (110) 

where -P(C) is the probability for the entropy created between — T and T to be ( given 
the protocol X(t) and the non-conservative force f([^],t). -P r (C) is the same probability, 
given the time-reversed protocol A and time-reversed force f r . 

Similar results can be obtained for isolated systems by switching off the interaction 
with the bath, i.e. by taking 7 = 0. It is also straightforward to obtain extended 
relations when the bath is taken to be out of equilibrium, for example by using 
r(t — t') 7^ j(t—t')+j(t' — t), and the contribution of the change in the dissipative action 
is taken into account. This kind of fluctuation relation may be specially important in 
quantum systems. 

5.2. Generic relations between correlations and linear responses 

A number of generic relations between linear responses and the averages of other 
observables have been derived for different types of stochastic dynamics: Langevin with 
additive white noise |30| . Ising variables with Glauber updates [31], or the heat-bath 
algorithm [321 1331 l3il 135] . and even molecular dynamics of hard spheres or Lennard- Jones 
particle systems [36] . Especially interesting are those in which the relation is established 
with functions of correlations computed with the unperturbed dynamics [301 E3] as 
explained in [37]. The main aim of the studies in [31] -[37] was to give the most efficient 
computational method to obtain the linear response in the theoretical limit of no applied 
field. Another set of recent articles discusses very similar relations with the goal of 
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giving a thermodynamic interpretation to the various terms contributing to the linear 
response [38]- [30]. 

In the concrete case of Langevin processes this kind of relations can be very simply 
derived by multiplying the equation by the field or the noise and averaging over the noise 
in the way done in [30] ■ We derive here the same relations within the MSRJD formalism, 
using a symmetry property that is more likely to admit an extension to systems with 
quantum fluctuations. 



5.2.1. A symmetry of the MSRJD generating functional valid also out of equilibrium. 
We consider the most generic out of equilibrium situation. We allow for any initial 
preparation (P;) and any evolution of the system (F). jT>[ip,ijj] e s ^'^ is invariant 
under the involutary field transformation T e0 m, given by 



T 

' c.a 






2/3 



r 



-i 



Eq„[ 



[111) 



M'(VO J v u - v M'(^) ' 

The meaning of the subscript referring to 'equation of motion' will become clear in the 
following. For additive noise [M'(ip) = 1] the transformation becomes 



\$ u 1-4 -iip u + 2/3 / T 



mip v 



+ / 'Jv-wlpv 



and in the white noise limit it simplifies to 



tyu !->> -hp u + /?7o X 



mijj u - F u [ip] +7o-0„ 



(112) 



The proof of invariance is similar to the one developed in Sect. I4.2l when dealing with the 
equilibrium symmetry. The Jacobian of this transformation is unity since its associated 
matrix is block triangular with ones on the diagonal. The integration domain of ip is 
unchanged while the one of ip can be chosen to be the real axis by a simple complex 
analysis argument. In the following lines we show that the action S evaluated in the 
transformed fields remains identical to the action evaluated in the original fields. We 
give the proof using an additive noise but the generalization to a multiplicative noise is 
straightforward. We start from the expression (TTUj) and evaluate 



S[T e0 mlp, Teomlp] = In P;(^_ T , ^_ T ) + 



it[j u - 2/3 / r-± w EQ 



x 



EQ M [^] - I J fl-^u-v (-^ w + 2f3JT- 1 _ z EQ z [ij} 



lnPi(V>_T,^-T) + 



r4 - 2/3 / T-\Eq v 



- I p 1 r u _ u , itp u 



(113) 



30 



Contrary to the equilibrium transformation 7~ eq , it does not include a time- reversal and 
is not defined in the Newtonian limit (7 = 0). 



5.2.2. Supersymmetric version. In Sect. 14.101 we encoded the fields ip, iip, c and c* in 
a unique superfield \I/. In this fashion, the transformation % m given in eq. f illip acts 

as 



2/3 f IT 1 u M'(y(u,9,9*))EQ u m 



d t M(S>(t,9,9*)) 
and leaves the equilibrium action S[$f], see eq. (1731) . invariant. 



(114) 



5.2.3. Out of equilibrium relations. We first derive some relations in the additive case 
\M'(ip) = 1] and then we generalize the results to the case of multiplicative noise. 

Additive noise. Using T e0 m m the expression (1351) of the self response R(t,t') we find 
{i/) t ^t')s = (T eom iptTeom.iipt>}s = -(i/>A)s + 2/3/ r~^ v (tp t EQ v [ilj]) s , 

J v 

giving an explicit formula for computing the linear response without perturbing field: 
R{t,t') =/3 UvT-\t' -v) 



x 



md 2 v C(t, v)+ du 7(1; - u)d u C{t, u) - (iP(t)F{[iP],v)) 



(115) 



Once multiplied by V t n_ t i and integrated over tl yields 
md?,C(t, t') + J du 7 (t' - u)d u C{t, u) 

- (ip{t)F{[ip],t')) = /3" 1 / du T(t' - u)R(t, u) , (116) 

with no assumption on the initial Pi(ip_T,4>-T)- 
Using now % om in (EQ t [ip}iip tl ) s , we get 

{Ed t [i)]iipt') s = (^Qt[%omtp]%o m ii>t'}s 

- (EQ t [V]i^)s + 2/3 fr^ u (EQ t [ij}EQ u [ij}) s . 

J u 

Since (Eq^jEq^ [ip]) s = /3 _1 T t _ u , this simplifies in 

(EQ t [V>]iVv)s = 5 t -t> , 
that yields 

md?R(t, t') + [ dv <y(t - v)d v R{v, t') - (if(lf)F{[i/j\, t)) s = S(t - t')(117) 
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with no assumption on the initial Pi. One can trade the last term in the lhs for 
P J u Ttf_ u {£(u)Ft[iJ)])(: by use of Novikov's theorem. 

Integrating both eqs. (11161) and fll 17j) around t — t' we find the equal-time conditions 

md t/ C(t,t%, =t = 0, m d t R{t,t%,^ t - = 1, m d t R{t,t%,_^ t+ = . (118) 

The last two relations imply that the first derivative of the response function is 



discontinuous at equal timegjj; 



The use of this symmetry is an easy way to get a generalization of eq. (11151) for a 
generic response Rab- Indeed, applying this transformation to expression (1361) of the 
linear response we obtain 



- d:{A[ m\Fm,u)f^)s 

+ Jdv 7 (u - v^HmtMv ^^ s} . (119) 

This formula gives the linear response as an explicit function of multiple-time correlators 
of the field ip. For example, if B is a function of the field only (and not of its time- 
derivatives), just the n = 0-term subsists in the above sum: 

R AB (t,t') = 0Jdu T-\t> - u) im %(A[mMu) d * [ f.u,h s 



d V(t') 



{mt)]Fm , n) a ^m )s 



+ yd,, 7 (« - v)a,(AW)Wv) a g [ *W ] h} ■ (120) 

As another example if one is interested in the self-response of the velocity, v4[-?/'(t)] = 
B[ip(t)} = d t ip(t), one obtains 

R AB (t,t') = (3 JduT-\t'-u) {m d t d 3 u C(t,u) - d t d u (4>(t)F([4>],u)) s 

+ fdv-f(u-v)d 2 v C(t,v)} . (121) 

Multiplicative noise. Similar results can be obtained for multiplicative noise. Applying 
the transformation to the correlator J T t > - u (iptM' (ifj t >) M' (ip u )itfj u ) $ we obtain 



<VtEQt#]>s = /r 1 / r t ,- u y> t M'(fo)M'u> v )u/> u )s , 

J u 

ft It is clear from the expressions given in (|118[) that the overdamped m — > limit allows for a sudden 
discontinuity of the response function as well as a finite slope of the correlation function at equal times. 
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implying 

md 2 C(t,t') + [ w-u{>l>tM , U> v )M , U> u )8 u i/> u )s 

J u 

- {IhFtrWhs = r 1 f Y t ^ u {^ t M'{^ tl )M'^ u y4 u ) s . (122) 

J u 

Applying now the transformation to the correlator (EQ t [if)]iip t ') s , one obtains 

(EQ t [V]i#)s = St-tf + P' 1 f V^ u (M\i) t )M'{^ u y4 u ^ v ) s (123) 

J u 

yielding 

md 2 R(t,f) + f lt _ u (M'^ t/ )M'(i/j u )d u ^ t/ ) s 

J u 

- (F t [iP]iip t ,}s = 5t-t>f3- 1 [ r t -. u (ip t M'(ip t ,)M'(ip u )i$ u ) s . (124) 



a 



One can check from eqs. (I122p and (I124p that the equal-time conditions given in 
eqs. (J118P are still valid in the multiplicative case. 



5. 3. Composition of % ora and %, 



cq 



For an equilibrium situation, the MSRJD action functional is fully invariant under the 
composition of % ora and T eq , 

that simply reads in the white noise limit 

For simplicity we only show the implication of this symmetry in this limit and in the 
additive noise case: 

R(t,t')= - R(-t,-t') + ^ [md 2 C(-t,-t') + (^-t)V'^(-t'))) s ] . 

7o 

Using equilibrium properties, i.e. time-translational invariance of all observables and 

time-reversal symmetry of two-time correlation functions of the field ip (shown in 

Sect. I4.2p . and causality of the response, 

R(r) = 0(t)£ [md 2 C(r) + <lK*)VMO)>s] , (125) 

7o 

with t = t — t' which is eq. ( I116p after cancellation of the lhs with the last term in the 
rhs when FDT between R and C holds [also eq. (11171) after a similar simplification]. 
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Here again, one can easily obtain a generalization of this last relation for a generic 
response Rab by plugging the transformation into the expression (I36l) of the linear 
response. 

6. Conclusions 

In this paper we recalled the path-integral approach to classical stochastic dynamics 
with generic multiplicative colored noise. The action has three terms: a deterministic 
(Newtonian dynamics) contribution, a dissipative part and a Jacobian. We identified 
a number of symmetries of the generating functional when the sources are set to zero. 
The invariance of the action is achieved by the three terms independently. 

One of these symmetries applies only when equilibrium dynamics are assumed. 
Equilibrium dynamics are ensured whenever the system is prepared with equilibrium 
initial conditions at temperature /3 _1 (a statistical mixture given by the Gibbs- 
Boltzmann measure), evolves with the corresponding time- independent conservative 
forces, and is in contact with an equilibrium bath at the same temperature /3 _1 . The 
invariance also holds in the limit in which the contact with the bath is suppressed, i.e. 
under deterministic (Newtonian) dynamics, but the initial condition is still taken from 
the Gibbs-Boltzmann measure. This symmetry yields all possible model-independent 
fluctuation-dissipation theorems as well as stationarity and Onsager reciprocal relations. 
When the field-transformation is applied to driven problems, the symmetry no longer 
holds, but it gives rise to different kinds of fluctuation theorems. 

We identified another more general symmetry that applies to equilibrium and 
out of equilibrium set-ups. It holds for any kind of initial conditions - they can be 
any statistical mixture or even deterministic, and the evolution can be dictated by 
time- dependent and/or non-conservative forces as long as the system is coupled to an 
equilibrium bath. The symmetry implies exact dynamic equations that couple generic 
correlations and linear responses. These equations are model-dependent in the sense 
that they depend explicitly on the applied forces. They are the starting point to 
derive Schwinger-Dyson-type approximations and close them on two-time observables. 
Although the symmetry is ill-defined in the Newtonian limit, the dynamic relations it 
yields can nevertheless be evaluated in the Newtonian case. 

Finally, we gave a supersymmetric expression of the path-integral for problems 
with multiplicative colored noise and conservative forces. We expressed all the previous 
symmetries in terms of superfield transformations and we discussed the relationship 
between supersymmetry and other symmetries. 

We intended to present a self-contained presentation of some symmetry properties 
of classical deterministic and stochastic dynamics. We focused on the so-called model 
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A dynamics (with a non-conserved order parameter) for a O-dimensional field. The 
generalization to a vector field is straightforward. Extensions to this work include the 
study of other dynamics such as the so-called model B dynamics (with a conserved 
order parameter). Higher dimensional fields would also allow the study of forces that 
do not work, such as Larmor precession around a magnetic vector field [6] . This article 
should serve as an introduction to and motivation for the study of quantum problems 
that we shall develop in [42J. Although some of these results were known, notably 
those associated to additive noise processes, they were scattered and somehow hidden 
in different publications. The close relation between all these properties was not always 
fully appreciated either. The multiplicative noise results are, as far as we know, new. 
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Appendix A. Conventions and notations 

is the Heaviside step function. When dealing with Markov Langevin equations, the 
choice of the value of the Heaviside step function 0(£) at t — is imposed by the choice 
of the Ito [0(0) = 0] or the Stratonovich convention [O(0) = 1/2]. However, away from 
the Markov case, i.e. as long as both inertia and the color of the bath are not neglected 
simultaneously, the choice of O(0) is unconstrained and the physics should not depend 
on it. We recall the identities 

f°° At f y 

/ — e ™v = S(y) and / dx 5(x) = G{y) , (A.l) 

J -oo ^ J —oo 

where 5 is the Dirac delta function. 



Field theory notations. Let ip be a real field. The integration over this field is denoted 
fT>[ip] . If A is a functional of the field, we denote it A[ip}. If it also depends on one or 
several external parameters, such as the time t and a protocol A, we denote it -AQ^], A, t). 
Whenever A is a local functional of the field at time t (i.e. a function of ip(t) and 
its first time-derivatives), we use the short-hand notation A[^(i)]. The time-reversed 
field constructed from ip is denoted ip: %p(t) = ip(—t). The time-reversed functional 
constructed from A([ip],X,t) is called A T : A r ([ip], X,t) = A([t/j], X, —t). Applied on local 
observables of ip, it has the effect of changing the sign of all odd time-derivatives in the 
expression of A. 
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To shorten expressions, we adopt a notation in which the arguments of the fields 
appear as subindices, ip t = if){t), It-t 1 = lit — t'), and so on and so forth, and the 
integrals over time as expressed as f. = J dt . 

Grassmann numbers. Let Q\ and #2 be two anticommuting Grassmann numbers and 
9\ and Q\ their respective Grassmann conjugates. We adopt the following convention 
for the complex conjugate of a product of Grassmann numbers: (6162)* = d^i- 

Appendix B. Discrete MSRJD for additive noise 

In this Appendix we discuss the MSRJD action for processes with additive colored noise. 

Appendix B.l. Discrete Langevin equation 

The Langevin equation is a stochastic differential equation and one can give a rigorous 
meaning to it by specifying a particular discretization scheme. 

Let us divide the time interval [— T, T] into N + 1 infinitesimal slices of width 
e = 2T/(N + 1). The discretized times are t k = — T + ke with k — 0, ..., TV + 1. The 
discretized version of ip{t) is ip k = ip(t k )- The continuum limit is achieved by sending 
N to infinity and keeping (N + l)e = IT constant. Given some initial conditions ip\ 
and ipi, we set ip\ = ipi and ipo = ipi — eipi meaning that the first two times (t and t±) 
are reserved for the integration over the initial conditions whereas the N following ones 
correspond to the stochastic dynamics given by the discretized Langevin equation: 

EQ fe _! = m 5 F k (ip k ,ip k -i,...) + e } j kl 

1=1 
= & , (B.l) 

defined for k — 1, ...,N. The force F k typically depends on the state ip k but can have 
a memory kernel {i.e. it can depend on previous states ip k -i, ip k -2, etc.). The notation 
7m stands for 7^ = e~ x f -du j(t k — ti + u). The £ k are independent Gaussian random 
variables with variance (£&&) = fi~ l T k i where Y k i = 1m + lik- Inspecting the equation 
above, we notice that the value of ip k depends on the realization of the previous noise 
realization ^ k _ 1 and there is no need to specify £ an d 6v+i- 

In the white noise limit, one has 7^ = e~ 1 i S k i, (£&&) = 2j f3~ 1 e~ 1 5 k i where 5 is 
the Kronecker delta, and 

EQ fc _! = m = F k (ip k , V'fe-i, ...) + 7o = 6c • 

e z e 
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Appendix B.2. Construction of the MSRJD action 

The probability density P for a complete field history (if; , ^ 1( ..., V'w+i) is se t by the 
relation 

P(V> , V'ij •••> V'J\r+i)dV'odV'i-d'0jv+i 

= P i (ip u iP i )dip i d->Pi ^n(6,6,-^iv)d6d6-d^ ■ 

P; is the initial probability distribution of the field. The probability for a given noise 
history to occur between times t\ and tjsr is given by 

Pniti,-,**) = M$ e-^M=^^ r M 6 (B.2) 

where T^ 1 is the inverse matrix of T^ (and not the discretized version of the inverse 

operator of T) and the normalization is given by M. 2 N = , / , _ n where det (...) stands 

aet fei(,P r fci J 

for the matrix determinant. From eq. ( 1B.2J) . one derives 

P(ij) 0l ip h ...,ip N+ i) = \JN\Pi{4>l, ^^)P n (EQ , ..., Eq^.O , (B.3) 

with the Jacobian 

j = det / 9fei^iv,^) \ = det / gjjMfe i EQo illll EQjv_i ) 
\d(ip ,7p 1 ,...,7p N+1 ) ) l I 9(^0,^l,---,^iV+l) 



that will be discussed in Appendix B.3 The expression (1B.2[) for the noise history 
probability reads, after a Hubbard-Stratonovich transformation that introduces the 
auxiliary variables ipk (k — 1, ..., A/"), 

A/WP n (£i,-,6v) = /d^ 1 ...d^ iV e- £ ^ i ^ fe+ 5 /3 " le2 ^ i ^ r ^ 

d^o-d^jv+i 5(^o)5(^+i)e- eEfci ^ EQfc - 1+ ^" l£2Efeii ^ r '=' i ^ , (B.4) 



with A/jv = (2' K / e ) N ■ in the last step, we replaced £& by EQ fe _! and we allowed for 
summations over k = and fc = A/" + 1 as well as integrations over ipo and tpN+i at the 
cost of introducing delta functions. The Hubbard-Stratonovich transformation allows 
for some freedom in the choice of the sign in front of i^ i n the exponent (indeed P n is 
real so P n = P*). Together with eq. ( IB. 31) this gives 

A/'ivP(?/'o,V'i,...,Vjv+i) = |JaH di/) Q ...dif) N+ i 8(i(>o)6(ii>N+i) 



xe -E fc i^EQfc-i + 2/ 3 " 1 E H i^r fci iV. i +lnP i (^ 1 ,^ tt ) 

that in the continuum limit becomes 

jVpUl _ Ij-r^ll glnPi / -prj,i g-Jdu mK«)Eq([</>],u)+± Jdudu i^(«)/3- 1 r(«-«)i'0(«) 



37 



with the boundary conditions ip(—T) = ip(T) = and where all the integrals over time 
run from — T to T. In the following, unless otherwise stated, we shall simply denote 
them by J . The infinite prefactor M = lim (2ir/e) can be absorbed in the definition 

N— s-oo 

of the measure: 

N N+1 

Vty,j>] = lim (i- ) TT<ty fc d& . (B.5) 

JV->-oo Kill/ - LJ - 



fc=0 



Markov case. In the Markov limit, the Langevin equation is a first order differential 
equation, therefore only the first time to should be reserved for integrating over the 
initial conditions. Moreover, one has to specify the discretization: 

tpk ~ tpk- 



EQ/c-i 



7o- 



- F k (ip k ) = £ k 



(B.6) 



where ip k = d^>k + (^~ o)ipk-i with a e [0,1]. a = corresponds to the Ito interpretation 
whereas a = 1/2 corresponds to the Stratonovich one (see the discussion in Sect. 12. 4p . 
Following the steps in Appendix B.2 we upgrade eq. ( IB. 61) to the following a-dependent 
actioqjj: 

tpk ~ ipk- 



S N (a) 



J2 l^io^kf-^k 



7o- 



F k (ip k 



- -KM 

7o 



(B.7) 



The last term in the RHS comes from the Jacobian: 

/<9EQfc_l\ = -j-r /To 



Jn = det 



A7 



o^(^* 



A' 



e -<E fc ^(^Jk) 



In the Ito discretization scheme (a = 0) this Jacobian term disappears from the action. 
Although Sn(o) seems to be a-dependent, we now prove that all discretization schemes 
yield the same physics by showing that the difference Sn(o) — SV(0) is negligible. The 
Taylor expansion of F k (i]j k ) around ^-l, F k (ip k _i) + a (ip k - ip k _ x ) F'(i{) k _ 1 )+0(e) [since 
^) k - ^ fc _i = 0(y/e)], yields 

S N {a)-S N {0) = aeJ2F'{ip k L ' 



>k-l, 



i4>k (ipk - i/>k-i) - 



7o 



+ 0(e 2 ). (B.8) 



Although the first term within the square brackets looks smaller than the second one, 
they are actually both 0(1) since iip k = 0{l/\/e). Thus, each term in the sum in the 
RHS is 0(e). We now compute the average of Sn(cl) — SW(0) with respect to S'jv(O) 
by neglecting in the latter the term e\il) k F k (%p k _i) which is of order y/e whereas the 
others are of order 1. Since (iip k (ipk — i>k-i)) Sif(o) = 1/7oj it is easy to show that 
(Sn(cl) — Sn(0))s n (o) = and therefore all the Sjv(o) actions are equivalent to the 
simpler Ito one. 



f We omit the initial measure which is not relevant in this discussion. 
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Appendix B.3. 



Jacobian 



Appendix B.3.1. Discrete evaluation of the Jacobian. In this section we take 

the continuum limit of the Jacobian defined in eq. f ]B.4j) . In the additive noise case, we 
start from 

0(^i,^i,EQ o ,...,EQjv_j 



Jn = det 



d{ip ,^i,---,^N- 



det 



-det 





dEQo 
dipo 

9EQ! 

dipo 

dipo 
9Eq 
djb 2 

•9EQ! 

dip 2 
9Eq^_! 



0... 

^EQi 
dips 



1 0... 

1/e 0... 

dEQ dEQ 

dipi dtp2 

dEQx dEQx 9EQ! 

<9i/>i dip2 



0... 

9EQ! 
9-03 



\ 





9Eq n _ 1 
dipN+i ) 



\ 



\ 



8ip2 





dip N+1 / 



(B.9) 



Causality manifests itself in the lower triangular structure of the last matrix. One can 
evaluate the last determinant by plugging eq. (IB.lj) . It yields 

(9EQ fc _! _ 1 (m\ N 



Ji 



N 



1 N 

in 



e JJ »fa« 



The Jacobian JT = lim J7iv is therefore a field-independent positive constant that can 
be absorbed in a redefinition of the measure: 



T>[ip, ip] = lim 



Af->oo e \27re 



m 



A? 



iV+1 



] J &ip k dtp k 



(B.10) 



fc=0 



We show that this result also holds for multiplicative noise in Appendix C 



Appendix B.3.2. Continuous evaluation of the Jacobian. One might also wish 

to check this result in the continuous notations. A very similar approach can be found 
in EHj. In the continuous notations, lim J N reads up to some constant factor 



J[i/j] = det, 



5Eq([7p],u) 
8ij){v) 



where det [...] stands for the functional determinant. Defining F^ v as SF u [if;]/Sip v , the 
Jacobian reads 



J[ip] = det, 



md 2 J> u - v + 



lu- 



, d,J. 



w K ~ / w u w—v 



F' 
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det t 
det t 

det, 



J w 
J w 



det, 



Ou—V I ^U—W-T will 

J w 



exp Tr u „ In [5 W _<, - M u 



exp - ^P - J \ M0M0...0M I (B.ll) 

n=l \ n times ) ,,,„ 



where we used the notations M uv = {G o F'} uv = J G tt _ w F^[^]. G is the retarded 
Green function solution to 



md^G(u — v) + / dw 7(14 — w)d w G(w — v ) = <5(w — v) . 



(B.12) 



Since both (?„_„ and F' uv are causal, it is easy to see that the n > 2 terms do not 
contribute to the sum in eq. ( IB. 111) . If the force F([ip],t) does not have any local term 
(involving the value of ip or ip at time t) the n — 1 term is also zero. Otherwise the 
n = 1 term can still be proven to be zero provided that G(t = 0) = 0. This will be true, 
as we shall show in the next paragraph, unless the white noise limit is taken together 
with the Smoluchowski limit (m = 0). Away from this Markov limit we establish 



J[ip] = det t 



mdl5(u - v) + / j u - w d w 8,. 



-W L/ tO ,J «) — V 



meaning that the Jacobian is a constant that does not depend on the field ip. 

We now give a proof that G(t = 0) = 0. Taking the Fourier transform of eq. ( IB. 121) 



G(t = 0) 



do; 

27 



G(u) = - 



du 



2-k mu 2 + iu^w) 



(B.13) 



G{uj) and 7(0;) are the Fourier transforms of the retarded Green function and friction. 
They are both analytic in the upper half plane (uhp) thanks to their causality structure. 
The convergence of the integrals around |a>| — > 00 in eq. (IB. 131) is ensured by either 
the presence of inertia or the colored noise. For a white noise [7(0;) = 70], it is 
clear that the mass term renders the integrals in eq. (1B.13J) well defined. In the 
m = limit the convergence is still guaranteed as long as the white noise limit is 
not taken simultaneously. Indeed, because 7(0;) is analytic in the UHP, it is hence either 
divergent on the boundaries of the UHP or constant everywhere [7(0;) = 70]. In the 
first case, which corresponds to a generic colored noise, this renders the integrals in 
eq. ( IB. 131) well defined. In the second case, corresponding to a white noise limit, they 
are ill-defined and require a more careful treatmentqj. When the integrals in eq. ( 1B.13I) 
are well defined on the boundaries, the absence of poles (or branch cuts) in the UHP of 



f In the white noise limit, G(t) = 7^ [l — e _7o */ m ] 9(t) is a continuous function that vanishes at 
t = 0. If we take m — > in the previous expression, we still have G(0) = and G(t) = 0(£)/7q for 
t 3> JTi/70. By choosing 0(0) = 0, these two results can be collected in G(i) = 0(i)/7o for all t. The 
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G(oS) gives, after a little deformation of the integration contour in eq. (IB. 13|) above the 
cu = pole, the result G(t = 0) = 0. 

Appendix B.3.3. Representation in terms of a fermionic field integral. The 

determinant can be represented as a Gaussian integration over Grassmannian conjugate 
fields c and c*. This formulation is a key ingredient to the supersymmetric representation 
of the MSRJD path integral. Let us first recall the discretized expression of the Jacobian 
obtained in eq. (IB. 91) : 

1 fdEQk^' 
Jn — -aetfci | - 7r> , , 



where k and I run from 1 to N. Introducing ghosts, it can be put in the form 
Jn = -^ /dc 2 dc*...dc JV+1 dc^_ 1 e e ^ fc =° U ^ c ^ °*i q 

- --^ /dc dc*...dc J v + idc* Ar4 . 1 cW^ce* -'-" ^'-" ; " ' ; 



where in the last step, we allowed integration over cq, ci, c;jv and c^- +1 at the cost 
of introducing delta functions (remember that for a Grassmann number c, the delta 
function is achieved by c itself). In the continuum limit, absorbing the prefactor into a 
redefinition of the measure, 

1 1 JV+l N+l 

VtyJ] = lim -—- jj- ] J dip k dj k and 2>[c, c*] = lim ] \ dc k dc* k , (B.14) 

TV— >oo [ZTT) £ J - J - TV— >oo - LJ - 

v y fc=0 fc=0 



this yields 

,7^]= fv[c,c*]e sJ ^ c *^ 
with 

J u J V 



' C"U 1 



b~i>v 

and the extra boundary conditions: c(—T) = c(—T) = c*(T) = c*(T) = 0. Plugging the 
Langevin equation (JTJ, we have 

6E Qu M_ _6FM r 



dip v Oip v 



■w 



Jacobian is still a constant. This limiting procedure where inertia has been sent to zero after the white 
noise limit was taken, is the so-called Ito convention. However if m is set to from the beginning, in the 
so-called Stratonovich convention with 6(0) = 1/2, then G(t) = Q(t)/jo for all t and G(0) = 1/(270). 
This can lead to a so-called Jacobian extra-term in the action. If F([i/j],t) is a function of ip(t) only 
(ultra-local functional), it reads — 1/(270) / F^(jp u ). It is invariant under time-reversal of the field 
Tpu l— > ip—u as long as F' is itself time-reversal invariant. 
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The kinetic term in S J [c, c*, ip\ can be re- written 

c* dl5 u - v c v = / c* d 2 u c u + ©o [c*c - c*c\^ T + 6 <5 [c*cf- T ■ 

The last two terms in the RHS vanish by use of the boundary conditions (c_^ = c_y = 
c* T = t* T = 0). The retarded friction can be re- written 

C* d u J u - v C v -Q C* [j u+ T C-T - lu-T Ct] , 
J u 

where the second term vanishes identically for two reasons: the boundary condition 
(c_t = 0) kills the first part and the causality of the friction kernel (7„ = OVm < 0) 
suppresses the second one. If there is a Dirac contribution to 7 centered at u = like in 
the white noise case, the other boundary condition {c*_ T = 0) cancels the second part. 
Finally, we have 

S J [c, c\ ^ = / < d 2 u c u + [ fcl \d ulu . v - *^M] c v . (B.15) 

Ju Ju Jv L OiPv 

Appendix C. Discrete MSRJD for multiplicative noise 

The discretized Langevin equation reads: 

^Q k -i = m = F k (ip k , ip k _ u ...) 

e z 

1=1 

with ipk = aipk + (1 — a)ipk-i an d k = 1,...,N. In the Markov limit (m = and 
7h = €~ 1 j 5m) the results depend on a (see the discussion in Sect. I2.4p . In the additive 
noise case, the choices a = and a = 1/2 correspond to the Ito and Stratonovich 
conventions, respectively. However, we decide to stay out of the Markov limit: the 
results are then independent of a and we choose to work with a — 1. The probability 
for a field history is 

P^ ,1pl,...,1p N+ l) = \J N \Pi(lpl, - ~ — )Pn(EQ , ..-, VQN-l) , (CI) 

where we introduced the shorthand notation EQ fc = Eq k /M'(ifjk+i)- The Jacobian is 

j _ dct f d(lpi,1pi,£l,---,£N) \ __ dct / g^j^i.EQo^.^EQjy.i) . 



P n is still given by expression ( 1B.4J) and P n (EQ , ..., Eqat-i) reads, after the substitution 
A^ 1 /d^-d^jv+i <K$o)*(Av+i) l^le-^^^^-^^-^Swi^w^jrwM'^M ; 
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where J7jv — det/y (4;M'(^)) is the Jacobian of the previous substitution. The 
probability for a given history is therefore 



JnJi 



N 



xe -E fc i4EQ fc _i + i/3~ 1 E H i^M'(^ fe )r fc; M'(V'i)i^+lnPi(V'i, :: ^ Q ) 

The Jacobian Jn defined in eq. (IC.2I) reads 

Jn - - det H (^^ ^- - ^^ EQ,,, <W 

" 7 ^ detfc/ ^^+7 ~ m 7 ^ EQ *- 1 Wl J ( } 

where fc and I run from 1 to N. Causality is responsible for the triangular structure of 
the matrix involved in the last expression. The second term within the square brackets 
yields matrix elements below the main diagonal and these do not contribute to the 
Jacobian. Therefore, we find 

N 



JnJn — - y | 



dEQ^ 1 /mN v 

fc=i 



e l= x 9ipk+i e Ve 2 . 

that is the same field-independent positive constant as in the additive noise case that 
can be dropped in the measure, see eq. (IB.lOj) . 

A fermionic functional representation of the Jacobian can be obtained by 
introducing ghosts, expression (IC.3I) can be put in the form 



with 



11/" e 

JnJn = — jj / dc dc* ...dc N+ idc* N+l c* N+1 c* N cic e' 



iV+liV+l -. _ N+l -mthi , \ 

fc=0 1=0 Yl k=0 KYK+LJ 



In the continuum limit it becomes 

S^ hm S*~ f [ c * 6Eq ^h - fc* M "^ Eo U\c 
S -^J N -JJ v Cu S^ ° v J u Cu M^ u ) EQu[ ^ C - 

with the boundary conditions c(—T) = c(—T) = and c*(T) = c*(T) = and the 
measure of the corresponding path integral is given in (1B.14[) . 
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